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Abstract. Ultra-high-energy neutrinos and cosmic rays produce short radio flashes through the Askaryan effect when they 
impact on the Moon. Earthbound radio telescopes can search the Lunar surface for these signals. A new generation of low- 
frequency, digital radio arrays, spearheaded by LOFAR, will allow for searches with unprecedented sensitivity. In the first 
stage of the NuMoon project, low-frequency observations were carried out with the Westerbork Synthesis Radio Telescope, 
leading to the most stringent limit on the cosmic neutrino flux above 10 23 eV. With LOFAR we will be able to reach a 
sensitivity of over an order of magnitude better and to decrease the threshold energy. 

Keywords: neutrinos, cosmic rays, radio telescopes, Moon, radio emission 
PACS: 95.85.Ry, 95.55 Jz, 95.55.Vj, 96.20.-n 



1. INTRODUCTION 

A century after their discovery, cosmic rays (CRs) still 
represent one of the unsolved questions in astrophysics. 
While modern experiments like the Pierre Auger Obser- 
vatory (Auger) [1], HiRes [2] and the Telescope Array 
(TA) [3] can efficiently detect ultra-high-energy (UHE) 
CRs, two complications remain that make it hard to iden- 
tify their sources. First, the trajectories of CRs are bent in 
the (inter-)Galactic magnetic fields, so the arrival direc- 
tions of the CRs at Earth do not align with the positions 
of their sources. This deflection depends on the charge of 
the CRs, and therefore on their composition. An Auger 
study of possible correlation with AGNs [1] disfavors CR 
isotropy above 6- 10 19 eV. However, composition stud- 
ies by Auger show a trend towards a heavier compo- 
sition at the highest energies [4], which suggests large 
separations between CR arrival directions and sources. 
In contrast, HiRes data does not reject isotropy and sug- 
gests a light CR composition [5]. A second complication 
is that CRs above an energy of 4 • 10 19 eV will interact 
with the cosmic microwave background (CMB). In these 
Greisen-Zatsepin-Kuzmin (GZK) interactions pions are 
created, which produce neutrinos when they decay [6, 7]. 
A steepening of the CR spectrum at this energy has been 
found experimentally [8, 2]. The attenuation length of 
protons due to the GZK effect is <~ 50 Mpc. Sources out- 
side this radius cannot be found by studying UHE CRs 



directly. The detection of UHE neutrinos that are pro- 
duced in GZK interactions or directly in the source is an 
attractive alternative to find these sources. Neutrinos are 
not affected by magnetic fields, and they can travel over 
cosmic distances almost unattenuated. 

It was first suggested by Dagkesamanskii and 
Zheleznykh [9] to use the Moon as a detector for 
UHE neutrinos and CRs. At energies above several PeV 
the Moon is opaque to neutrinos. Neutrinos impinging 
on the Moon will interact below the lunar surface. In 
a charged current interaction <~ 20% of the energy is 
converted into a hadronic cascade. The leptons that are 
created are in principle detectable, though in the case of 
muons and taus the relevant cross-sections at UHE are 
highly uncertain. Electrons initiate an electromagnetic 
shower which is hard to detect with the lunar Cherenkov 
technique because it is elongated due to the LPM effect 
[10]. However, at sufficiently high energy, the length and 
depth of the cascade is reduced because of the increase 
of the photonuclear and electronuclear cross sections 
[11]. Particle cascades beneath the lunar surface can be 
detected by the radiation they emit through the Askaryan 
effect [12]. 

GLUE [13] was one of the pioneering experiments at 
high frequencies (>1 GHz). Presently, LUNASKA [14] 
performs high frequency measurements with ATCA and 
the Parkes radio telescope. At the Expanded Very Large 
Array, RESUN [15] performs 1.4 GHz lunar observa- 



tions. 

The NuMoon project observes in a low frequency win- 
dow which offers an optimal detection probability at the 
highest energies. In the first phase, observations were 
performed with the Westerbork Synthesis Radio Tele- 
scope (WSRT) leading to the most constraining limit on 
the UHE neutrino flux above 10 23 eV [16, 17]. Presently, 
we are preparing observations with the Low Frequency 
Array (LOFAR) [24]. 

2. RADIO EMISSION FROM LUNAR 
PARTICLE CASCADES 

Neutrinos and CRs impacting on the Moon initiate cas- 
cades below the surface, which develop a time-varying 
negative charge excess propagating at the speed of light 
in vacuum, and thereby produce radio emission [12]. Al- 
though the radiation is most intense at the Cherenkov an- 
gle, it can be misleading to think of it in terms of pure 
Cherenkov radiation. The rise and decay of a moving 
charge itself would also produce radiation if the index of 
refraction of the medium were unity [18, 19]. Actually, 
at low frequencies, when the wavelength is larger than 
the cascade length, the radiation is coherent at angles far 
from the Cherenkov angle [20] . 

The spread of the radiation around the Cherenkov 
angle is of large importance for the detectability of the 
radio pulse. To observe the pulse on Earth, it needs to 
be able to escape the Moon. Since the Cherenkov angle 
is equal to the angle of total internal reflection on the 
surface, the high frequency component can only leave 
the Moon if the shower points upwards to the surface. 
This is only the case when a particle enters the Moon 
from behind, near the lunar rim. The low frequency 
component, on the other hand, is spread out over a much 
larger opening angle, and part of the radiation can escape 
the Moon even if the cascade is directed downwards into 
the Moon. At ^150 MHz the whole visible surface of the 
Moon contributes to the effective detector volume [21]. 

We simulated the distribution of observable pulses 
over the lunar surface for an isotropic neutrino flux with 
an energy spectrum E~ l in the range 10 21 — 10 24 eV. We 
assume that 20% of the neutrino energy is deposited in a 
hadronic cascade for all flavors. The neutrino cross sec- 
tions [22] are extrapolated to high energies for a 500 m 
deep layer of regolith, and the radio pulse parametriza- 
tion of Alvarez -Muniz et al. [23] is used for the hadronic 
showers. As a mean value for the attenuation length for 
the radiated power we have taken X — (9/v[GHz]) m. 
Surface roughness is not included since it does not play 
an important role at low frequencies [21]. 

The geometry definitions used here are sketched in 
Fig. 1: 9 V is defined as the angle between the direction of 



V 




FIGURE 1. Geometry: 6 V is defined as the angle between 
the direction of the neutrino and the observer; b is the relative 
distance of the point where the pulse emerges from the lunar 
surface to the center of the face of the Moon. 
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FIGURE 2. Distribution of radio pulses with power greater 
than 500 Jy as a function of energy and b. 
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FIGURE 3. Distribution of radio pulses with power greater 
than 500 Jy as a function of 8 V and b. 



the neutrino and the observer; b is the relative distanceof 
the point where the pulse emerges from the lunar surface 
to the center of the face of the Moon (b = 1 corresponds 
to the rim of the Moon). 

In Fig. 2 the distribution of radio pulses above a 
threshold of 500 Jy in the 1 10- 190 MHz band is plotted 
as a function of neutrino energy and b. The energy spec- 
trum has been reweighted to E~ 2 and the total observed 
flux has been normalized to unity. The color scale thus 
represents the relative probabilities of observing an event 
in different parts of the parameter space. The threshold 
energy for detection is ~ 10 21 eV and at one decade in 
energy higher the whole surface of the Moon contributes 
to the detector volume. Even at higher energies events 
are more likely to occur close to the rim. This is mainly 
a projection effect: the high-Z? region represents a very 
large part of the visible half of the lunar sphere. 

In Fig. 3 the same set of events is now plotted as a 
function of 9 V and b. There is a strong clustering of 
events towards the rim of the Moon with an arrival angle 
9 V in the range 110° - 150°. Neutrinos from 50° - 90° 
can produce observable pulses on the complete lunar 
surface, while for even smaller values of V , observable 
events are rare. Hence, localization of a low-frequency 
radio pulse on the Moon provides constraints on the 
arrival direction of the neutrino. 



3. NUMOON WITH LOFAR 

LOFAR [24] is a new kind of radio telescope consisting 
of thousands of simple omni-directional antennas that 
takes advantage of recent developments in fast electron- 
ics. The received radio signal is digitized at each antenna 
and sent to a central processor where directed beams are 
synthesized by applying the appropriate time delays. Be- 
cause the antennas themselves contain no moving com- 
ponents they are relatively cheap and can be constructed 
in large amounts. Another great advantage of synthesiz- 
ing beams through software is that multiple beams in dif- 
ferent directions can be in operation at the same time. 

Each LOFAR station consists of 96 Low Band An- 
tennas (LBA), which operate in the 10—80 MHz range, 
and 48 High Band Antennas (HBA), which cover the 
1 10—240 MHz band. Currently, the array consists of 24 
core stations, concentrated in an area of ~12 km 2 , 9 re- 
mote stations, distributed over the northern part of the 
Netherlands with baselines up to 80 km from the core, 
and 8 international stations. 

For the NuMoon experiment we will use the HBA an- 
tennas operating in the 110— 190 MHz window. The 
beamforming is done in two stages because of limited 
bandwidth between stations and the central processor. 
At each station a beam of a few degrees opening an- 
gle is formed in the direction of the Moon. This sig- 




FIGURE 4. Arial view of the LOFAR core. The island, 
called the superterp contains 6 of the 24 core stations. It is 
divided in larger and smaller patches which house the LBA and 
HBA antennas respectively. 



nal is transported to the central processor where it is 
first corrected for ionospheric dispersion. At LOFAR fre- 
quencies, a typical ionospheric electron content of 10 
TEC-units (TECU=10 16 electron/m 2 ) causes a disper- 
sion which spreads a bandwidth-limited signal pulse over 
hundreds of nanoseconds [25]. In order to be able to per- 
form a pulse search the signal needs to be de-dispersed 
with an accuracy of 1 TECU. An interesting approach is 
measuring the absolute TEC value via the Faraday ro- 
tation of polarized light due to ionospheric plasma and 
the Earth magnetic field. We are currently investigating 
the possibility to use the polarized light of the rim of the 
Moon for this purpose. 

In the next step, the de-dispersed station beams are 
combined to form array beams of ^0.1 degree width. 
Up to 50 such beams can be synthesized simultaneously 
to cover the complete surface of the Moon. Each of these 
beams will be searched in real-time for short pulses. 
When a candidate pulse is found in one of the beams a 
coincidence check will be performed. A radio pulse orig- 
inating from a certain location on the Moon is strong 
in only one of the array beams (although neighboring 
beams may also receive a signal, depending on beam 
shape and sidelobes). This is a powerful vetoing tech- 
nique that removes pulses from strong noise sources near 
one of the stations: a very strong noise pulse in one of the 
station beams will cause a spike in all the array beams. 

Each of the LOFAR antennas is connected to a ring 
buffer that stores the last five seconds of raw data. When 
a pulse is found at the central processor that passes the 
anti-coincidence veto, a trigger will be sent back to the 
stations and the buffers are read out and permanently 
stored. This allows us to use the full resolution, full 
bandwidth raw time-series data of all antennas for offline 
analysis [26]. 



4. SENSITIVITY 



5. CONCLUSIONS 



In the past years, much progress had been made in im- 
proving the sensitivity to cosmic neutrinos over a large 
range of energy, with a large range of techniques. In 
Fig. 5 existing upper limits on the all-flavor neutrino flux 
are plotted in solid lines. The most stringent limits are 
set by IceCube [27], ANITA [28] and the WSRT [16] at 
increasingly larger energies. Also plotted are limits pub- 
lished by Rice [29] and Auger [30]. 

The expected sensitivity that is reached with one week 
of LOFAR data will improve the previous WSRT limit by 
an order of magnitude, while also decreasing the thresh- 
old detection energy. In the future, the Square Kilometer 
Array (SKA) [31] will provide an even better sensitiv- 
ity using the same techniques as LOFAR. At lower en- 
ergies, better sensitivities will be reached by future ex- 
periments ARA [32] and ARIANNA [33] (not plotted), 
both searching for radio flashes of neutrinos in ice, and 
the JEM-EUSO satellite [34]. 

The flux of GZK neutrinos is sensitive to various un- 
known parameters including source evolution, CR com- 
position and source energy cutoff. Figure 5 contains pre- 
dictions computed by [35] based on varying input param- 
eters. The shaded band represents a region which con- 
tains a large part of the possible models. The model that 
produces the largest neutrino flux is indicated by dash- 
dotted line. 

UHE neutrinos can also be created in the decay of 
supermassive particles (top-down scenarios). The TD- 
curve in Fig. 5 delimits the parameter space in which 
neutrinos may be found originating from moduli — 
weakly coupled scalar particles predicted in supersym- 
metric theories — from the kinks of cosmic string loops 
[36, 37]. 

A cascade induced by a CR will be initiated just below 
the lunar surface in contrast to neutrino induced show- 
ers. Therefore, the radio signal suffers less attenuation. 
However, the question rises whether or not the radiative 
field can properly develop in a volume that is smaller 
than the wavelength of the radiation. A detailed deriva- 
tion [18, 19] shows that such a formation zone is not 
needed. This becomes evident when one realizes that a 
system of growing and decaying current would also ra- 
diate if it were to develop in a vacuum. Hence, the di- 
electric medium is not essential for the formation of the 
radiative field, although it certainly affects it. In Fig. 6 
the WSRT upper limit on CR flux [18] is plotted, as well 
as expected sensitivities for LOFAR and SKA (low band 
and middle band). The data points represent the high end 
of the CR spectrum as measured by Auger [8]. 



The Numoon experiment uses the Moon as a ultra-high- 
energy cosmic -ray and neutrino detector. It is different 
from other lunar Askaryan experiments as it searches for 
radio flashes at lower frequencies (1 10 — 190 MHz). This 
has the advantage that a much higher sensitivity can be 
reached, although the threshold energy is slightly lower 
at high frequencies (compare for example the sensitivi- 
ties of the SKA low and middle band in Fig. 6) . Observa- 
tions with the WSRT have led to the most stringent upper 
limits above 10 23 eV. Preparations are now in progress to 
perform measurements with LOFAR. The expected sen- 
sitivity is over an order of magnitude better than previous 
measurements. In the future, the SKA will make another 
big leap in sensitivity using the same techniques. 
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FIGURE 6. CR flux limit of WSRT, as well as expected 
sensitivities of LOFAR and SKA. The dotted line indicates 
a constant particle flux (cm _2 sr _1 s _1 ) corresponding to the 
particle flux of the highest-energy data point. A flux above this 
line would have been found by Auger. 



